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FACILITY: MATH LIBRARY 
ABSTRACT: 


This module contains the routine MTHSHMOD: 
It returns the roneineer of the division of argi/arg2 using 


the age e 
. Tee eri a sare2d ears? 


AUTHOR: Jeffrey C. Wiener, CREATION DATE: 21-DEC-1982 
MODIFIED BY: 
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O38 -SBTTL HISTORY ; Detailed Current Edit History 
009 ; 3-00 Original version of complete re-write JCW 21-DEC-82 
4 : DIVIDEND changed to equal 8 and DIVISOR changed to equal 4. JCW 14-Jun-83 


nm 


"E=SE=1984 91:35:93 EATARTE. saeamrunmopemar:1 P29 By 
.SBTTL DECLARATIONS 

INCLUDE FILES: 
NONE 

EXTERNAL SYMBOLS: 


DECLARATIONS 


Sete Se Ge Ge Ge Ge 


-DSABL GBL :; Force all external symbols to be declared 
~EXTRN MTHSSSIGNA 
-EXTRN MTHSK *FLOUNDMAT 
~EXTRN MTHSK_ INVARGMAT 
; LIBRARY MACROS CALLS: 

SSFDEF ; Define SF$ (stack frame) symbols 
EQUATED SYMBOLS: 


P.112 = *x00000070 ; 112*2°0 
IGA MASK = FOIFF 


OWN STORAGE: 
NONE 
PSECT DECLARATIONS: 
.PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 
; CONSTANTS: 


TWO_EXP_112: 
LONG “*X00004071, “XO, “XO, “XO ; 2e*112 
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Suet -SBTTL MTHSHMOD - H REAL*16 remainder 
: FUNCTIONAL DESCRIPTION: 
| 
| 
| 


Return the remainder of He H aaa point format 


Remainder = argi - (int(argl/arg2))*arg 


The algorithm used to evaluate the HMOD function is as follows: 


X = the first argument. 
Y = the second argument. 
step 1. m = the exponent of Y. 
n = the exponent of X. 
c=n-™m 
If c < 0, end with result = X. 
step 2. 1 = the fractional part of X. 
J = the fractional part of Y. 
If lor J, lele-d 
Go to ae P 
step 3. L = 2“(p-1)*I, where p = 113 for H_floating numbers. 
step 4. T = ae 
| = Bs delat T is int(L/J) or int(L/J)+1 
If 1<0, 1leI1¢4J T was int(L/J)+1 
step 5.c =c = (p-1) 
If ¢ > 0 go to step 3. 
step 6. If c = oa go to step 9. 
step 7. L = 2*(p-ltc) * I 
step 8. l=L-J#T 
step 9. Result = 2m * | 


CALLING SEQUENCE: 

CALL MTHSHMOD (remainder.wh.r, dividend.rh.r, divisor.rh.r) 
INPUT PARAMETERS: 

The two input parameters are H_floating-point values. 


Bete Ge Se Se Ge Se Be Ge Ge Se Ge Se Se Ge Se Ge Se Se Se Ge Ge Se Ge Ge Se Se Ge Ge Ge Ge Se Ge Se Ge Se 


DIVIDEND = : Dividend = X in the algorithm. 

DIVISOR = 12 ; Divisor = Y in the algorithm. 
IMPLICIT INPUTS: 

NONE 


OUTPUT PARAMETERS: 


Remainder is the remainder of the division of 
argi/arg2, returned as an H_Floating point value. 


IMPLICIT OUTPUTS: 
NONE 
COMPLETION CODES: 
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54 4000 8F 
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7E 54 6€E 


ot 
MTHSHMOD - H REAL*16 remainder “SE 


— 


SIDE EFFECTS: 


+ 


16909 09 09 C909 0D 09 G9. 09 SINISE NSN DPA PAA AA AA IIIT 
WOONAUE WN O OO NAU EWN 0 OD NAMES WN OC OONOAULS WN OUOan 


SOSOSSSSSSSSSSSOSSSSSSSSSSSSSOSSSSSSSSSSOSSSSSSSOSO 
BRR RR RRR RRR RRR R EE RWWA 2 ae 


D NNN F NNW WN AANA FT OH OOW MM AMM ONINVIOOOOOOOOOCOOo 


Ooooeo 


POCO 
ce ee ee ee ee ee ed ed ed od od ed dd 8 = 8d wd 8 wd od od 2 — 2 2 3 3 2) — 2 — — — 3 I 


SOCOoOOOOOWOOOOO 


SOOOOOOOOOOSOOSOOOSOOSSOSOOSOSOSOSOSOOOSOOSOSOSOSOOOSOVOOSOOOOOOOOOOOOOOOOOOoOO 


OOO fLOOUO 
WS OCOnNOuUlwn—o 


rRorornore 


i=) 
< 
m 
z 


P=1984 01:38:0 AX/VMS Macro V04-00 P 4 
Pa} 9b4 91:38:99 EMTHRTL. SRCOATHHMOD. MAR; 1 ee (3) 


NONE 


MTHS_INVARGMAT = Invalid argument to math Library if the divisor is zero. 
MTH$_FLOUNDMAT - Ae ag underflow in math Library is signaled if 
the FU bit is set in the callers PSL. 


-ENTRY MTHSHMOD, “M<R2, R3, R4, RS, Rb, R7, RB> 

MOVH Q@DIVISOR(AP), R4 3; RG/R7 = Y 

BEQL ERROR ; Y=0. Division by zero 

MOVH Q@DIVIDEND(AP), RO 3; RO/R3 = X 

BICL #*XFFFF8000, R4, -(SP) ; (SP)=m is the biased exponent of Y 
BICL #*XFFFF8000, RO, RB ; R8B=n is the biased exponent of X 
SUBL2 (SP), R8 ; R8 = ¢ = n-m unbiased 

BGEQ TEP ; Result is X if X<Y, ie, if c<0 
MOVO RO, 34(AP) 3 @4(AP) = X 

RET 

MOVW #*xX4000, R4 3 R4/R7 = J = biased ifract(Y): 
MOVW #*x4000, RO ; RO/R3 = I = biased ifract(X): 


In STEP_4 and STEP_8 the calculation of I = L = J*int(L/J) must be 
computed as precisely as possible. To do this we will need to write J as 


Jed 
where J1 = the high 56 bits of J and J2 = J - J1, the low 57 bits of J. 
Bi Oy naae is used to extract the 7 bits of J from longword3 that belong 
) 


MOVO R4, -(SP) 3 (SP) = J 

BICL  #HIGH_MASK, 8(SP) : 

CLRL 12(SPJ ; (SP) = J1 replaces the value 
: of J on the top of SP 

SUBH3 (SP), R4&, -(SP) 3; (SP) = J2eJ- J 

CMPH RO, R4 s 1f I<J 

BGEQ 1 Gea y 

BRW STEP_5 i go to STEP_S 

SUBH2 R4, RO ; else I = I-J 

BLEQ OVER 

BRW STEP_5S ; go to STEP_S if 1>0, or 


else the algorithm ends 
TSTW @DIVIDEND (AP) ; the sign of the result is 


MTHSHMOD 16-SEP-1984 01:38: AX/VMS Macro V04-00 P 
3-002 MTHSHMOD - H REAL*16 remainder aioe 0:38:09 EMTHRIL. SREIMTHHMOD MAR: 1 ~— (3) 
4 C 4 BGEQ DONE 3; the same as the sign of 
§ 72FD 5 MNEGH RO, ; the first argument, A. 
04 BC 7DFD f § DONE: peye RO, a4(AP) ; Return the correct result 
7 
50 01 . ig 5 ERROR: ASHQ #15, #1, RO ; Y=0. Reserved operand 
2 7 7C 10 CLRQ R 
7E OO'8F QA 7 11 MOVZBL #MTHSK_INVARGMAT, (SP) 3; error code 
00000000 ' GF 9 fe 4 19 CALLS #1, G*ATHSSSIGNAL ; signal the error 
04 BC 0 ror? 4: 1? peye RO, @4(AP) ; Return the correct result 
§ 15 
50 00000070 8F CO St 18 STEP_3: ADDL2 #EXP_112, RO 3 RO/R3 = L = 2e2(p-1)*] 
0096 18 
B38 5 
0096 1 ; STEP_4: 
0096 ¢ 3 2*(p=1) = 2°(112) is added and then subtracted from 
Opoe Z 3 T = int(L/J) to ensure that T = chopped(L/J) or chopped(L/J)+1 
os ae 
7E 54 7DFD 0096 5 MOVO R4, -(SP) ; save_J for use in STEP_5S 
54 50 54 67FD OO9A § 8 DIVH3 R4, RO, RG 3: R4/R7 = T = L/ 
54 FFSC CF 60FD O09F 9 ADDH2 TWO_EXP_112, R4 > RG/R7 = T = Tease(pot) 
54 FF56 CF 62FD BOAB $39 SUBH2 TWO_EXP_112, R4 j aay ty L/J chogeed 
3 or choppe 
Q0nB 232 “ 
Sonn 534 5 
00AB 3¢ : The calculation of I =L =- J*int(L/J) must be computed as precisely 
ae § 3 as posstole. Je dé this we will need to write T as 
a 8 : where 21 = the high 56 bits of T and 72 = T = 21, the low 57 bits of T. 
ear rt : Now, using J = J1 + J2, 
00AB 4g F L-J * int(L/u) = Lb = (1 + U2) ® (21 + 22) 
QOAB 43; =LtL-=- (21 * J1) = (21 * 43} 
OAB 44; - (22 * J1) _ (22 * J2) 
OAB 45 ; =tL-(21*5)2 (22% J) 
QOAB $6 3 
yt 
7E 54 7DFD AB 49 MOvVO R4, -(SP) 3; Stack Z = INT(L/J) 
O8 AE FFFFOIFF BF CA OOAF 250 BICL  #HIGH_MASK, 8(SP) : Start to form 21 
AE D4 7 51 CLRL 1 ; (SP) = 71 
7E 54 6€ O37 A 26 SUBH (SP), R4&, -(SP) : (SP) = 22 
54 10 AE 40 AE 65FD F 5 MULH 64 (SP), 16(SP), RG 3 Comp te Z1*J1 
4 62rd 90ce 334 SUBH + RO 3 = ts gteut 
54 30 AE 10 AE 65FD OOCA 5 MULHS 16¢SP), 48(SP), R4 > RG/R7 = 118) 
50 54 oghd D1 2$ SUBH 3 Rove st- qteJ 
54 40 A E 65FD 0005 MULHS (SP), 64(SP), R4 > RG/R? = 12%) 
5 ogre DB 8 SUBH : Rone zsLe- pied = 22*J1 
54 30 AE 6 65FD OODF MULHS = (SP), 48(SP), R4 ; R4/R7 = Z2*J 
5 4 62FD O0ES 260 SUBH2 4, RO : RO/R3 = L = ZeJ 
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af 0 co —9 61 ADDL2 #32, SP remo e a 72,J1,J2 from stack 
4 7DED i 6¢ nov (SP5+, R4 t ROR 
; : 4 BGTR STEPS ; If RO/R3>0 the algoritm cotinues 
19 F4 5 BLSS SUBTRACT_J 
069 31 4 § BRw RETURN ; If RO/R3=0 the algorithm ends 
oF 67 SUBTRACT_J | 
50 54 60FD Fg $8 ADDH2 = R4, RO ; Add J back in because you had 
are 3 ; T = chopped(L/J)+1 
58 00000070 B9 i6 orp 4 STEP_S: suai 2 eexe -j12. R8 ie = co(p-1) = c-112 | 
58 00000070 C 166 rs ADDL2 #EXP~112, RB 3 ¢ = (pel)+c = 112¢4¢ 
B88 Bre i 
3198 7 : The next two Lines of code are STEP_6 and STEP_7. 
4 
67 13 010D Hi BEQL STEP_9 : 
50 58 CO 010F B6 ADDL2 R8, RO 3 Lb = [*2*(c+t) 
He e $ 
011 85 ; 2*(p-1) = 2*(55) is added and then subtracted from | 
011 86 ; int(L/J) to ensure that T = chopped(L/J) or chopped(L/J)+1 | 
S113 5a | 
o BR HE Boe Bm ae 
54 FEEO CF 60FD 0118 291 ADDH Fu ni At, R4 > R4/R7 rst = Te2eepe1) 
54 FEDA CF 62FD 0121 3 SUBH TWO-EXP_112, R4 : ity -}) =LL/Jc opped 
0127 9 s /J chopped + 
0127 94 ;+ 
0127 95 ; 
3 4 38 3 STEP_8: 
0127 298: The calculation of I = L = J#int(L/J) must be {computed as precisely 
3 + : as posers rhs do this we will need to write T as | 
a 4 4 : where 71 = the high 56 bits of T and 22 = T = 21, the low 57 bits of T. | 
3 4 § : Now, using J = J1 + J2, 
127 5: L-J * int(L/u) = Lb = (il + U2) ® (21 + 22) 
127 BS 3 =L = (21 * J1) = (21 * J2) | 
127 07 ; - (22 * JI). te: * J2) 
127 4 : =t-(21"J5)2+(22%/)) 
127 : 
eg 
7E rey 7DFD 127 \¢ MOVO R4, ~(SP) 3 Stack Z = INT(L/J) 
08 AE ererone? +4 ca : R \7 art wHIGn MASK, 8(SP) ; rier te, fore 
7E 56 «O6E 63FD 136 315 SUBH3 (SP), R4, -(SP) : (SP) = 22 | 
54 10 AE 40 AE 65FD 01 18 MULHS © 64(SP), 16(SP), R4 : Compute 21#J1 | 
50 54 62FD 1 1 SUBH R4, RO 3; RO/ =L = 21*J1 


54 SO AE 10 i 65FD 
50 5 ? FD 

iad ae 
54 30 A . o3fe 

i Bae 

me 

6 14 

1A 13 

50 54 60FD 


20 AE 4000 BF A 
50020 13 ‘ 


B5 
18 
50. 8000 8F AB 

04 BC 0 vere 


04 8 
0D 04 AD 06 


Q00000000°8F  odD 
00000000'GF 01 fe 
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16(SP), 48(SP), R4 


@DIVIDEND(AP) 
RETURN 
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a4 (AP) 
#SFSV_FU, SFSW_SAVE_PSW(FP), 


#MTHSK FLOUNDMAT 
#1, G*ATHSSSIGNAL 


pelelelelelelelelelelelelelelelelelelalelelelelol ele elel el elo a) oa) 
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Sete Ge Ge Se Se See 


Sete Se Be 


NO_Fu 


(ee ] 


R4/R7 = 21*J2 
RO/R3 =L = Zi*J 
= 


R4/R tJ 
RO/R3 =L- fied = 22*J1 

R4/R7 = 22*J 

RO/R3 - I*J 


~ 
Remove 1 and 72 from the stack | 
Restore J | 


If RO/R3=0 the algorithm ends 
Add J back in because you had 
T = chopped(L/J)+1 


Remove the bias from m 
and form RO/R1 = 2*m*L 


the sign of the result is 
the same as the sign o 
the first argument, X. 
Return the correct result 


Set up default result to 0.0 


Branch if caller has not enabled F 
Report MTH$_FLOUNDMAT 

Signal the condition 

Return 


The working it Limit was 900 pages. 

7106 bytes ‘ 4 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space of located to hold 51 non-local and 0 local symbols. 
51 source has were read in Pass 1, producing 14 object records in Pass 2. 

pages of virtual memory were used to define macros. 


16- 1: AX/VMS Macro V04-00 P 
g-sEF 71982 te 38: 9 ci THRTL.SR SMTHHMOD . MAR; 1 7 ) 
MTHSHMOD 
MTHSK_FLOUNDMAT teeeeeee x 
MTHSK_INVARGMAT eeeeeeee = § a 
NO_FU 
VER 
RETURN 0 
FSV_FU 
SFSW_SAVE_PSW = 
TEP_ 0 
He 0 
STEP_ 0 
TEP 9 0 
SUBTRACT_J 0 
TEST_SIGN 0 
TwO_ExP_112 0 
ERFL 0 
wm omnm cw em mama moecre + 
! ;_Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
- ABS . 94.4 44444 ‘ 00 ( 0.) NOPIC USR CON ABS ~ LCL NOSHR NOEXE NORD wont NOVEC BYTE 
SABSS 0000 0:) 01 ¢ 1.) NOPIC USR CON ABS LCL NOSHR EXE RD T NOVEC BYTE 
_MTHSCODE BO0001A9 42 02 ¢( 2.) PIC USR CON REL LCL SHR’ EXE RD NOURT NOVEC LONG 
enna 
' Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 9 00:0 00-18 ae iy 
pomnend processing ! 1 09: : gt 0:0 
Syabol iable sort 0 if 7 HB 8:38 4 
cess f table output 0: : 00; :0 
Psect synopsis output 3: 0: 3: 6:4 
Cross-reference output :00:00. Ba: 2-08 
Assembler run totals 354 :00: : 0:00:14.6 


1< 


HSHMOD I 5 
vAKet Macro Run Statistics "O78EB =}38¢ 9}: 38: 9 MTARTL sae SATHGPOD OMAR: 1 Pic dy i. 


tower ce mee wes emma wem een moana + | 
H Macro Library statistics H 
Macro Library name Macros defined 
-$255SDUA28: CSYSLIBISTARLET.MLB; 2 4 
88 GETS were required to define 4 macros. 
There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHMOD/OBJ=OBJ$:MTHHMOD MSRC$:MTHHMOD/UPDATE=(ENHS:MTHHMOD) 
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